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ABSTRACT 



Aims. The influence of a time- varying ram pressure on spiral galaxies in clusters is explored with a new simulation method based on 
the N-body SPH/tree code GADGET. 

Methods. We have adapted the code to describe the interaction of two different gas phases, the diffuse hot intracluster medium (ICM) 
and the denser and colder interstellar medium (ISM). Both the ICM and ISM components are introduced as SPH particles. As a galaxy 
arrives on a highly radial orbit from outskirts to cluster center, it crosses the ICM density peak and experiences a time-varying wind. 
Results. Depending on the duration and intensity of the ISM-ICM interaction, early and late type galaxies in galaxy clusters with 
either a large or small ICM distribution are found to show different stripping efficiencies, amounts of reaccretion of the extra-planar 
ISM, and final masses. We compare the numerical results with analytical approximations of different complexity and indicate the 
limits of the Gunn & Gott simple stripping formula. 

Conclusions. Our investigations emphasize the role of the galactic orbital history to the stripping amount. We discuss the contribution 
of ram pressure stripping to the origin of the ICM and its metallicity. We propose gas accumulations like tails, filaments, or ripples to 
be responsible for stripping in regions with low overall ICM occurrence. 

Key words. Galaxies: general - Galaxies: interactions - Galaxies: intergalactic medium ~ Galaxies: clusters: general ~ ISM: structure 



1. Introduction 

Galaxies in clusters and rich environments are observed to be 
stripped of their interstellar medium, which quenches subse- 
quent star formation. The stripping can be due either to tidal 
interactions, spiral galaxy mergers into ellipticals, or to ram- 
pressure stripping from the intracluster gas. Both tidal and ram- 
pressure stripping are always tightly linked. Gunn & Gott (1972) 
assume that, after the formation of a galaxy cluster, the remain- 
ing gaseous debris is thermalized via shock heating to virial tem- 
peratures corresponding to random motions in the cluster, e.g. to 
a few times lO^K. This hot plasma at densities of a few times 
10"^ cm"^ influences the ISM in disks of spiral galaxies and can 
remove part of their gas through the ram pressure induces by the 
galaxy motions through the ICM. 

The ram pressure stripping is difficult to model, since several 
complex gaseous phenomena are involved. The first aspect is the 
simple pressure force, and this can be modeled with a simple 
algorithm using ballistic and sticky particles to represent the gas 
(e.g. Vollmer et al. 2001). But the full hydrodynamical processes 
include thermal evaporation, turbulent and viscous stripping, and 
also outflows due to star formation. 

The first notice of the 3D N-body/SPH simulations with the 
gravity tree to mimic the dynamical effect of the ram pressure 
on galaxies in clusters is by Kundic et al. (1993). The same ap- 
proach has been adopted by Abadi et al. (1999). The ICM is 
represented as a flow of particles along a cylinder of radius 30 



kpc and thickness 10 kpc. A spiral galaxy is in face-on, edge-on, 
or an inclined orientation relative to flow of ICM particles. With 
the simulations they examine the radius up to which the ISM is 
removed and compare the results with the prediction of Gunn & 
Gott (1972) that the ISM is removed from the disk if flie ICM 
ram pressure exceeds the restoring force: 
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where picM is the ICM density, v the relative velocity of the 
galaxy and the ICM, <^(r, z) the total gravitational potential of 
the galaxy as a function of the galactocentric cylindrical coordi- 
nates, and ZisM the ISM surface density. The conclusion is that 
Eq. [T] applies in the face-on galaxy orientation when the bulge 
does not dominate the disk gravity. In other cases, edge-on ori- 
entation and in central parts of galaxies, where bulge dominates, 
the stripping is less efficient. The time scale for gas removal is 
~ 10** years. 

The ISM is far from continuous, but instead consists of sev- 
eral components: hot, warm, and cold medium. Hot ionized and 
warm neutral HI medium are much more continuous compared 
to cloudy, cold, and molecular component. The SPH description 
relates more to the continuous ISM components, while molec- 
ular clouds are not represented well in the SPH simulations. 
Vollmer et al. (2001) adopt another approach, using N-body sim- 
ulations with sticky particles representing the inelastic collisions 
between molecular clouds. The ram pressure is included with an 
additional friction force acting on clouds in the wind direction. 
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Vollmer et al. (2001, 2006) introduce a time dependent ram pres- 
sure, which corresponds to the variations of the ICM density and 
of relative ICM-ISM velocity along the galaxy orbit in the clus- 
ter. They allow a moving ICM, which may increase its velocity 
relative the ISM to values higher than 4000 km s ' . 

With an isothermal SPH gas model, Abadi et al. (1999) show 
that gas stripping can be quite efficient in the core of rich clus- 
ters on a time-scale of 10^ yrs. The most efficient is face-on ori- 
entation. With a Eulerian code, Quilis et al. (2000) probed the 
efficiency of viscous coupling, enhanced by the presence of HI 
deficiencies in the center of galaxies. The ram pressure stripping 
does not strongly depend on the vertical structure and thickness 
of the gas disk (Roediger & Hensler 2005). In more detail, the 
effect of inclination of the moving galaxy is important as long as 
the ram pressure is comparable to the gas pressure in the galaxy 
plane. In general the effect is similar for all inclinations except 
for edge-on (Roediger & Briiggen 2006). The orientation of the 
gas tail behind the galaxy is not a good tracer of the galaxy mo- 
tion on its orbit. Acreman et al. (2003) simulate the infall of an 
elliptical galaxy and its gas stripping by the ICM showing that 
the formation of an X-ray wake at the first passage should be 
observable. 

The principal competing mechanism for perturbing galax- 
ies in clusters and stripping their gas is tidal interactions. High- 
speed galaxy encounters in clusters, galaxy harassment events, 
are efficient mechanisms for disturbing galaxy disks, and pro- 
viding stars to the intracluster space (Moore et al. 1996). Also 
low-speed galaxy encounters in groups before the formation of 
the cluster are able to drag tidal stellar tails out and lead to galaxy 
mergers and formation of ellipticals, which are then known to be 
more numerous in clusters. The evolution of encounter debris in 
galaxy haloes has been followed by Mihos (2004). 

Intracluster light is a consequence of the tidal stripping of 
galaxies. It has been determined that intra-cluster stars are older 
than galaxy stellar populations and more centrally concentrated 
in the cluster. The fraction of stars in the ICM is increasing with 
the richness of the cluster and it is usually a few percent of all 
stars, up to 20% (Arnaboldi et al. 2003, Murante et al. 2004). 
If the ICM metallicity can be explained partly by gas strip- 
ping from galaxies (Domainko et al. 2006), it can also come 
from intra-cluster supernovae exploding in the ICM. The irregu- 
lar structures in galaxy clusters found in X-ray are tracing cold 
fronts and shocks due to galaxy interactions or to the infall of a 
galaxy group. 

The fate of the stripped ISM in clusters is not clear, because it 
can end very hot and become part of the ICM or stay quite cold 
condensing into clouds with a molecular core and an HI outer 
boundary. Oosterloo & van Gorkom (2005) observed a large 
HI cloud near the center of the Virgo cluster and suggest that 
this cloud has been ram-pressure stripped from the galaxy NGC 
4388. The dense clumps in the HI plume might be molecular, 
and star formation could occur precisely at these high-density 
places. 

In this paper, we focus on computations of the ram pressure 
stripping, to evaluate its efficiency more precisely and to follow 
the fate of the stripped gas. The particular goal is to take the fi- 
nite time for this stripping interaction into account, as the galaxy 
passes quickly through the central region of the cluster. Most 
computations until now were carried out in the hypothesis of a 
stationary wind of hot gas. But the ram pressure acts only in a 
short lapse of time, as simulated by Vollmer et al. (2001). After 
this short period, a certain fraction of the gas falls back onto 
the galaxy. These latter simulations did not take into account 
both the pressure forces and the hydrodynamical physics of the 



ISM-ICM interactions. Roediger, Briiggen & Hoeft (2006) have 
followed the wake of gas produced by a quasi-stationary wind; 
here, on the contrary, we emphasize the impulsive character of 
ram pressure stripping. 

Assuming a static ICM described with a /3-profile, we ex- 
plore the ISM-ICM interaction in detail during the galaxy cross- 
ing of the central part of the cluster. The mass-loss rates are ex- 
amined as a function of the galaxy type, size, and mass. We are 
able to follow the gas stripped from the galaxy, forming a gi- 
ant tail of material. We tackle the issue of the origin of the ICM 
and the actual role of gas stripping in building it as a function of 
the cluster mass and richness. The importance of gas stripping 
could account for the metal enrichment of the ICM as a function 
of cluster type (e.g. Domainko et al. 2006, Schindler et al. 2005). 
After presenting the model of stripped galaxies and of the cluster 
in sections 2 and 3, we introduce the initial conditions in section 
4 and describe our simulation method in section 5. In section 6, 
we discuss the results of the numerical simulations, the simula- 
tion test is shown in section 7, and the equation of motion of an 
individual ISM gas element and the impulse approximation are 
given and applied to stripping in section 8. The present simula- 
tions are compared to other simulations in section 9. Section 10 
presents the discussion and conclusion. 

2. The galaxy model 

Our model of a spiral galaxy is a standard three-component halo 
-I- bulge -I- disk configuration. The halo and bulge components 
are given with spherically symmetric, Plummer density distribu- 
tions 

^bh 

Phjiid) = po —: ^-TTT, (2) 

where po = ^M^jJAna^i^ is the central density, d the distance 
from the galactic center, Mh, M/, are total masses, and Ob and a/, 
are the radial scaling factors corresponding to the bulge or halo. 

The axially symmetric disk follows an infinitely thin 
Toomre-Kuzmin disk multiplied by a sech^(z/zo) term that de- 
fines its isothermal vertical profile, 

Pd(r,z) ^po . sech2(z/zo), (3) 

(r2 -t- fl2)3/2 

where (r, z) are the galactocentric cylindrical coordinates, and 04 
and zq are the disk scaling factors. The central density is deter- 
mined by Po = Md/Ana^zo with the total mass Mj of the disk. 

We introduce model spiral galaxies of various types: late- 
type (L) or early-type (E), subdivided according to their total 
mass: massive (M) or low-mass (m). The L-types have very 
low bulge-to-disk mass ratios compared to E-types, and M-types 
have a massive halo component compared to m-types. Thus we 
employ four models: LM, Lm, EM, and Em. The values of the 
model parameters are given in Table [1] The scale height of the 
disk in all models is zo - 0.25 kpc. 

Our standard galaxy model is the LM-type. Its rotation curve, 
the escape velocity, and the gravitational restoring force fresr = 
2^iSM X d<i>{r,z)/dz are shown in FiglT] where <l)(r, z) is the total 
gravitational potential including all components of the galaxy, 
and d<i>/dz the total force attracting ISM elements towards the 
galaxy's symmetry plane. The escape velocity is obtained as a 
function (r, z) from the potential: Vesc('", z) = -\J2\<^(r, z)\, and its 
maximum is always in the symmetry plane z = 0. Rotation ve- 
locities of other three models, Lm, EM, and Em, are about 160 
km s"', 220 km s"', and 150 km s"', respectively. 
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Fig. 1. LM model galaxy. Left: Rotation curve with contributions of the halo (dotted), bulge (dash-dotted), and disk (dashed) components. Center: 
Escape velocity v^^c as a function of the z-distance behind the disk plane for 16 radii r starting with r = i kpc (top curve) and a spacing of 1 kpc. 
Right: Gravitational restoring force /,,,.„ as a function of z for 16 values of radius r from 1 kpc (top curve) to 16 kpc. Maxima in z for individual 
radii are denoted with circles. 
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Fig. 2. Left: Maximum restoring force over z as a function of the galactocentric distance r for the four galaxy models adopted: LM (dashed), Lm 
(full), EM (dotted), and Em (dash-dotted). Center: Positions of the maximum values of the restoring force behind the disk plane for different radii 
r. Right: Surface density Sism of the ISM as a function of the galactocentric distance r. 
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Table 1. The disk, bulge, and halo parameters for late (L) or early (E), 
massive (M) or low-mass (m) type models of a spiral galaxy. 



mass. The ISM surface density Eism is plotted as a function of 
the galactocentric distance for the four galaxy types in the right 
hand panel of Fig.|2l 

3. The model of a galaxy cluster 

The hot ICM was observed with X-ray satellites ROSAT and 
ASCA and later resolved spatially with satellites XMM-Newton 
and Chandra (Bohringer 2004). The large-scale ICM distribu- 
tion in clusters is described with a j6-profile (Cavaliere & Fusco- 
Femiano 1976; Schindler et al. 1999): 



At a given r, the restoring force has a maximum at some 
distance z from the symmetry plane. With increasing radius r, 
this maximum is farther and farther from the z = plane and 
its values decrease. In Fig. |2] (left panel) we plot the maxima of 
the restoring force as a function of r for different galaxy types. 
The LM-type galaxy has a massive disk and massive halo, which 
gives the highest restoring force compared to other galaxy types. 
Both the disk and halo are less massive for the Lm-type, con- 
sequently, the restoring force is smaller; only in the central part, 
where the very concentrated bulge dominates, the restoring force 
increases rapidly. The bulge is much more massive and extended 
in an EM-type galaxy, so the restoring force in the outer galaxy 
parts overcomes that of an Lm-type. For radii between 1 and 4 
kpc, galaxy models Lm and EM show very similar runs of the 
maxima of the restoring force with r, but the maxima occur at 
lower z in the Lm-type (see Fig.|2] central panel). The Em-type 
galaxy, where the disk and halo have lower masses compared to 
the other types, has the weakest restoring force. 

We assume that the ISM initially follows the density dis- 
tribution of Eq. [3] and that it amounts to 10% of the total disk 



PiCM - PO.ICM 



p2 

'^c.ICM J 



(4) 



where picM denotes the ICM volume density at distance R from 
the cluster center, po.iCM is the volume density of the ICM in 
the cluster center, /?c,icm ^ parameter of the ICM central con- 
centration, and jSicM the slope parameter, which we set accord- 
ing to Schindler et al. (1999) to 1/2. Our standard cluster model 
has po.icM = 4 lO"-' cm"-' and ^c,icm = 13.4 kpc. Vollmer et 
al. (2001) use a ten-times higher value of po.icM to represent the 
ICM distribution in the Virgo cluster We discuss the difference 
in po.iCM setting in Sect. 19.31 

To model either a rich galaxy cluster with a lot of hot ICM 
or a poor cluster with just a little ICM in the center, we vary the 
values of ^c,icm and po.icM- We multiply either standard value 
with factors of 0.25, 0.5, 1, 2, and 4. The left panel of Fig. [3] 
illustrates the effect of varying Rc,vcm while po.icM stays fixed. 
The width of the ICM density peak grows with increasing /?c,icm- 

When the galaxy flies through the cluster, it encounters a 
certain amount of ICM. The right panel of Fig. [3] corresponds 
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Fig. 3. Left: Central view of the ICM density distribution in clusters 
with /?c.icM equal to 0.25, 0.5, 1, 2, 4 x 13.4 kpc, and po.icM = 4 lO"'' 
cm"'. Right: Peaks of the ICM density distribution with the following 
combinations of the j8-profile parameters (Rcicm, Poicm): (0.25, 3.82), 
(0.5, 1.95), (1, 1), (2, 0.53), and (4, 0.28) x (13.4 kpc, 4 lO^'cm-'). A 
galaxy traversing such clusters on a completely radial orbit from 1 Mpc 
distance crosses the same amount of the ICM. 

to five clusters with combinations of Rcjcm, Po.icm that make 
the galaxy cross the same amount of ICM when moving on a 
completely radial orbit, starting from a distance of 1 Mpc. The 
shape of the peak changes from narrow and high to wide and 
low. 

The gravitational potential of a galaxy cluster is produced 
by the distributions of ICM and dark matter (DM) in the cluster. 
The DM is introduced through another j6-profile. We set Rc,dm - 
320 kpc, po,DM = 1.6 10"^ cm"^, and /3dm - 1 as in Vollmer 
et al. (2001). The different steepness of gas and galaxy density 
profiles (J3icm versus /3dm) corresponds closely to observations, 
with the ICM density profile steeper than the galaxy distribution 
in the inner part of the subclump. This implies DM of either 3.7 
10'^ Mo or 1.4 10'"* Mq within 140 kpc (see later) or 1 Mpc from 
the cluster center Correspondingly, depending on the type of the 
cluster, ICM represents from 0.05% to 30% of the mass within 
the central 140 kpc (see TableO. 

4. Initial conditions 

The three components of the model galaxy are represented by 
a number of particles of equal mass: 12 000 for the halo, 6 000 
for the bulge, and 12 000 for the stellar disk. Their initial distri- 
butions follow the density profiles of Eqs.|2]and[3] The compo- 
nents are cut off at 40 kpc, 4 kpc, and 16 kpc, in the case of the 
halo, bulge, and disk, respectively. The ISM is represented with 
12 000 SPH particles of total mass equal to 10% of the total disk 
mass. This means that the mass of an ISM particle is in the range 
of (1.5 - 5.7) 10^ Mq in the four introduced galaxy types. Their 
initial distribution follows the density profile of the stellar disk 
(Eq.[3l see right panel of Fig.|2]). 

The ICM of the cluster is modeled with 120000 SPH par- 
ticles following the j6-profile (Eq. |4|i with a cut-off at R = 140 
kpc from the center of the cluster The total mass of the ICM 
particles is derived from cluster parameters. In the case of the 
standard cluster, there is 6 lO'" Mq of the ICM within the central 
140 kpc. Consequently the mass of an individual ICM particle is 
~ 5 10^ Mq, which is comparable to the mass of an ISM particle. 
For the LM-type galaxy and different clusters, the mass ratio of 
ISM to ICM particles ranges between 30 and 0.04. 

The cut-off radius of the ICM particles distribution set to 140 
kpc corresponds in the case of the standard model to the radius 
at which outer parts of the galaxy arriving from the cluster pe- 
riphery start to be affected by steeply rising ram pressure of the 




50 lOD 200 400 600 500 



Fig. 4. Radial profiles of the ICM density, pressure, and speed of sound 
about 0.5 Gyr before the galaxy enters the ICM distribution. The dashed 
curve in the density plot corresponds to Eq.|4] Lower right: Initial profile 
of the ICM velocity dispersion ctr, Eq.|5] 



crossed ICM. It means that we restrict the active sphere of the 
ICM distribution to central regions of the cluster, where the ef- 
fect of the ICM-ISM interaction is strongest. At larger distances 
only the gravitational effect of the ICM to the galaxy is consid- 
ered. 

A more realistic simulation would spread the ICM particles 
throughout the cluster so that the galaxy runs all its orbital time 
through the medium. However, as noted by Abadi el al. (1999), 
ICM particles that are too massive may become bullets punching 
holes into the ISM disk. To prevent this, the masses of individual 
ICM particles should at most be comparable to that of ISM par- 
ticles. With this assumption and the number of ISM particles set 
to 12000, the sphere of 1 Mpc radius of the cluster would need 
about two million particles of the ICM. Then the calculations 
would become too time-consuming and complex. 

Velocity distribution of the ICM particles is calculated from 
second moment of the collisionless Boltzmann equation provid- 
ing isotropy of the model (see e.g. Hernquist 1993): 

ai^— rp(R)^dR, (5) 
" p(R)Jr ' dR ' 

where 0,0, represents the total potential of the cluster's DM and 
ICM, and p {R) is the ICM density. The ctr profile is shown in 
the lower right panel of Fig. |4] Equation |5] is in fact an equa- 
tion for the hydrostatic equilibrium of the cluster After a slight 
initial relaxation, the system reaches a stable state, well before 
the galaxy enters the ICM distribution. Figure |4] depicts radial 
profiles of the ICM density (compared to the analytic profile of 
Eq. lU, pressure, and speed of sound about 0.5 Gyr before the 
galaxy arrives. To keep ICM particles within the 140 kpc radius, 
periodic boundary conditions are introduced such that, when a 
particle reaches the edge, it continues its course at centrally op- 
posite positions, with preserved velocity vector. Particles of the 
ICM are treated adiabatically, while particles of the ISM isother- 
mally with Tism = 10"^ K. 

Figure |5] shows the initial configuration of the system. 
Galaxy in face-on orientation relative to the orbit is located 
at the cluster periphery at 1 Mpc distance from the center, 
R - {x,y,z) - [0,0,1] Mpc, with no systematic velocity. 
Consequently, it falls freely towards the cluster center Thus, 



Jachym, Palous, Koppen & Combes: Gas stripping in galaxy clusters: a new SPH simulation approach 




-1000 -500 500 1000 960 980 1000 1020 1040 



Z [kpc] Z [kpc] 

Fig. 5. Initial configuration. Left: Galaxy falls face-on from 1 Mpc 
distance (dashed). The cut-off radius of the ICM is 140 kpc (dark), 
which outside (light) it operates only gravitationally. Right: Zoom on 
the galaxy with DM halo (dots) and disk+bulge (shaded) components. 
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Table 2. Parameters of the standard simulation run. Times at which 
galaxy enters the ICM and reaches the cluster center are stated, together 
with its velocity in the cluster center. 



in this work, all the orbits are completely radial, free-fall or- 
bits from the distance R = 1 Mpc to the cluster center For the 
standard cluster, the galaxy reaches the cluster center after about 
1.64 Gyr, at a velocity of about 1300 km s It means that it 
moves very supersonically through the ICM in the central re- 
gions (cf. Fig. m, and we expect a strong bow shock to form 
in front of the galactic disk. Domainko et al. (2006) draw at- 
tention to the role of the bow shock: the ram pressure behind 
the shock is lower since the relative velocities of the ICM parti- 
cles are smaller there than in front of the shock. On the contrary, 
Rasmussen et al. (2006) notice that, while the ram pressure is 
reduced at the shock, P,i, + pv^, where f is the thermal ICM 
pressure behind the shock, stays conserved for an inviscid fluid, 
and thus the reduction of the ram pressure behind the shock is 
balanced by an increase in static thermal pressure. Therefore, 
the force per unit area acting on the disk is close to the pre-shock 
ram pressure. 

Table |2] summarizes the parameters of the standard simula- 
tion run. Along the radial orbit, we can compute the ram pres- 
sure prain(R) = Picmv^,,; acting on the galaxy as it crosses the 
ICM. Ram pressure always peaks in the cluster center where the 
ICM density is highest, and also where the orbital velocity has 
its maximum. Following Fig. [3] for broader clusters the increase 
in the ram pressure starts at larger distances from the cluster cen- 
ter. With the same central density pojcM the peak value of p, „,„ is 
slightly higher for broader clusters due to the higher contribution 
of the ICM to the overall cluster gravity (Fig.[TTl left panel). 

In the simulations, we follow the galaxy on its orbit dur- 
ing 2 Gyr From the very beginning at 1 Mpc distance, the 
galaxy itself evolves and thus instabilities and spiral arms form. 
Consequently, the galaxy enters the ICM sphere at about 1 .5 Gyr 
as dynamically evolved. 

5. Simulation method 

We performed 3D tree/SPH N-body simulations using GADGET 
in version 1.1 (Springel et al. 2001). The tree method is a highly 
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Fig. 6. SPH sizes of the ISM (black) and ICM (grey) particles at the 
epoch when the galaxy runs through the cluster center. Mcm = 120 000, 
NisM = 12 000. Left: Face-on zoom view of the galaxy surroundings 
with the ISM particles located within ± 2 kpc about the disk plane, and 
the ICM particles within ± 20 kpc about the disk plane. Right: Edge- 
on view of the situation with ICM particles located within ± 20 kpc 
about the disk in the line of sight direction. Only every tenth particle is 
displayed. 

effective Lagrangian technique for computing gravity forces 
in non-homogeneous systems. The particles are aiTanged into 
groups hierarchically ordered from the smallest ones contain- 
ing only one particle up to the biggest one embracing the whole 
system. The scheme of an oct-tree structure dividing the compu- 
tational space into a sequence of cubes is used (Barnes & Hut 
1986). Then, the gravitational force exerted by distant groups on 
a particle is approximated by their lowest multipole moments, 
which reduces the computational costs to order 0(N log N). 

The smoothed particle hydrodynamics (SPH, Lucy 1977, 
Gingold & Monaghan 1977) is a Lagrangian method for solving 
hydrodynamical problems in which real fluids are represented 
by a set of fluid elements, i.e. particles. The key idea of the SPH 
is that every particle has an effective radius within which hydro- 
dynamical interactions with neighboring particles operate. Thus 
the particles are no longer point masses but are smoothed out, 
and local hydrodynamical properties are interpolated from the 
local mass distribution. The smoothing sizes of particles are de- 
termined to correspond to a radius within which a given number 
of neighboring particles is contained. 

GADGET in version 1.1 resolves only one gaseous phase, 
and all gas particles occurring in the system are hydrodynami- 
cally treated as one group. However, the volume densities of the 
ISM in the disk and of the suiTounding ICM differ by several 
orders of magnitude, as well as their temperatures. To treat their 
interaction correctly, the number densities of the ICM and ISM 
particles, nicM and hism, should be comparable. Otherwise, the 
ICM particles located far from the galaxy would be larger in the 
SPH sense than the ISM particles, since the number of neigh- 
boring particles is kept constant throughout the system. When 
these large ICM particles would later approach the disk, their 
sizes would be re-evaluated from the local mixing of ICM and 
ISM particles, and artificial punching effects could follow. 

To achieve similar number densities of both gases, a huge 
number of ICM particles are necessary. With these constraints in 
mind, we introduce a set of changes into GADGET related to the 
fact that we need to simulate two different gaseous phases (ICM 
and ISM) with different spatial resolution. First, the subroutine 
for searching the neighboring particles is adapted to distinguish 
between the two phases: the smoothing radii h's of either ICM or 
ISM particles are calculated separately from the neighbors of a 
corresponding phase. In other words, the search continues until 
Nh proper neighbors are found, while those of the second phase 
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Fig. 7. Time sequence - a galaxy crossing tiie cluster center. From top down: Edge-on projection and face-on cut (|z| < ± 2 kpc) of the ISM 
particles; the tail of the stripped ISM in the central distribution of the ICM particles. The time since the beginning of the free fall of the galaxy 
from the cluster edge is given in Myr in the first row. 



are omitted. Then a reasonable number of ICM particles, that 
are consequently substantially larger than the ISM, can be used, 
and their size does not change even when interacting with the 
disk. It ensures that all the disk area is covered by the ICM par- 
ticles. Second, the hydrodynamical interaction of the two phases 
can proceed: Both phases are searched together to find all par- 
ticles with \rjj\ < max{hi,hj). A given ISM particle then inter- 
acts with all the ICM particles within which smoothing radii it 
is contained. Third, the adiabatical and isothermal treatments of 
the ICM and ISM are implemented, respectively. 

Figure |6]compares the ISM (black circles) versus ICM (grey 
circles) particle sizes in the surroundings of the galaxy run- 
ning through the center of the standard cluster. Face-on and 
edge-on snapshots are displayed, and the A^icm = 120000 and 
NisM -12 000. Only every tenth particle is plotted, which gives 
the impression of voids occurring among particles. One notes 
that the size of the ICM particles is kept constant in the disk 
vicinity. 

To achieve a sufficient spatial resolution in the disk, we fix 
the number of ISM particles A^ism to an adequate value. In all 
simulations, the resolution in the disk is then the same. From 
Fig. |6]it follows that the ICM particles are substantially larger 
than their ISM counterparts, each of them covering a large part 
of the disk. That means that a group of close ISM particles lying 
within the smoothing radius of a given ICM particle feel only the 
ram pressure of the ICM but no pressure gradients, as would be 
the case if the sizes of ICM and ISM particles were comparable. 
In sectionl7l(Fig.[T8]l. we discuss test runs with varying numbers 
of ICM particles. 

Our simulation approach using only a reasonable number 
of ICM particles, and thus neglecting the detailed hydrodynam- 
ics in the ISM-ICM interaction treatment, has some similarities 



with the method of Vollmer et al. (2001), who in their sticky- 
particle simulations include the effect of the ram pressure only 
analytically as an additional acceleration on the ISM clouds lo- 
cated at the surface of the gas distribution. In our case, of course, 
the interaction between ICM and ISM is not restricted to the sur- 
face layer 

Our approach is more realistic than Vollmer et al. (2001) as 
far as the ICM is concerned. The ICM phase is treated fully hy- 
drodynamically. Therefore, as discussed later, simulations de- 
velop a bow shock that forms in the ICM in front of the galaxy. 
Consequently, the incoming ICM particles are deflected from 
their direction to flow around the disk. Moreover, when (from 
the SPH point of view) small ISM particles are stripped out of 
the disk, their number density falls and their size increases. In 
the tail of stripped particles, the sizes of both the ISM and ICM 
then become comparable, and the hydrodynamical treatment of 
the tail is correct (see Fig.|6]l. 



6. Simulation results 

First we describe results of the standard simulation run (see 
TablelU of the LM galaxy. This shows how the ISM of the galaxy 
changes under the influence of a) the ram pressure due to rela- 
tive motion of the ISM and ICM, b) the gravitational field of the 
cluster that changes along the orbit, and c) the evolving galactic 
gravitational field that includes the self-gravity of the ISM disk. 
We note that no eff'ects of the ram pressure to the stellar disk, 
bulge or halo galaxy components are observed in the simulations 
described below. 
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Fig. 8. The standard simulation run. Time evolution of the ISM mass 
within 16 different disk radii (1-16 kpc) and |z| < 1 kpc (dashed 
curves). The uppermost curve corresponds to the ISM mass bound to 
the galaxy. The instant of crossing the cluster center is marked with the 
vertical line. 



6.1. LM galaxy in the standard cluster 

Figure |7] displays snapshots of the ISM disk at five epochs of the 
galaxy's crossing through the central parts of the cluster Freely 
falling from the 1 Mpc distance, it approaches, at time 1 .52, the 
edge of the ICM particle distribution, where ram pressure rises. 
At 1 .64 Gyr it passes through the cluster center and experiences 
the peak ram pressure. At 1.76 Gyr it leaves the ICM on the 
other side of its distribution, relaxes after the fade-out of the ram 
pressure, and further continues towards the second apocenter. 
Figure |2] shows edge-on and face-on views of the ISM disk, and 
a wide view of the central parts of the cluster with the ICM par- 
ticles. 

When the galaxy enters the ICM distribution, the local ram 
pressure affects only the outer disk regions. The ISM from these 
regions is stripped and thus forms the end of the tail of stripped 
gas. Later on, more inner parts also get stripped by the increasing 
ram pressure, filling the central portion of the tail that is clearly 
visible in Fig.|2] The opening angle of the tail changes with time, 
which reflects the time evolution of the ram pressure vs. restor- 
ing force contest. Schulz & Struck (2001) describe one effect of 
the ram pressure enhancing the spiral structure of the disk since 
particles from inter-arm positions with lower surface densities 
are more easily removed. This effect is observed at time 1.58 
Gyr. 

The process of stripping is displayed in Fig. |8] It shows the 
number of ISM particles, i.e. the ISM mass, within different disk 
radii growing from 1 kpc to 16 kpc with a step of 1 kpc, and 
within a |z| < 1 kpc belt about the disk plane as a function of 
the orbital time. The uppermost curve shows the ISM mass Mbnd 
that is bound to the galaxy, i.e. the number of ISM particles with 
negative total energy. Before the galaxy enters the ICM zone 
at ~ 1.53 Gyr, no important changes are observed. Once the 
ram pressure rises, particles from decreasing radii are shifted 
out of the evaluation belt. With a delay of about 20 Myr after 
the galaxy's passage through the center, a maximum amount of 
the ISM is shifted out of the |z| < 1 kpc zone, and curves in 
Fig.[8]show a minimum corresponding to 49% of the initial total 
ISM mass. Particles down to about 4 kpc are influenced. After 
the minimum, a strong reaccretion of the ISM previously only 
shifted out of the disk by a short-term effect of the strongest ram 
pressure, but staying bound to the galaxy occurs. Of course, both 
the minimum value and the amount of the reaccreted material 
depend on the definition of the evaluation zone. 

With a constant flow of the ICM, Roediger & Hensler (2005) 
identify three phases of stripping: (i) an instantaneous phase 



when the outer disk parts are bent in the wind direction but stay 
bound to the galaxy; (ii) a dynamic intermediate phase, during 
which the bending breaks, a part of the ISM is stripped, and an- 
other part falls back to the disk; and (iii) a quasi-stable continu- 
ous viscous stripping phase when the outer disk layers are peeled 
off by the Kelvin-Helmholz instability at a rate of ~ 1 M© yr 
The duration of the first phase is ~ 20 - 200 Myr, while the 
following second phase is about ten times longer. In our simula- 
tions, we observe behavior similar to the first two phases; how- 
ever, the evolution of the Mbnd(0 and M|.|<ikpc(0 under the oper- 
ation of a peaked ram pressure is different (see Fig. [8]). Once the 
ram pressure rises, the outer disk layers are released, while the 
inner ones are only shifted from the |z| < 1 kpc zone. Although 
the rate of the release is lower compared to the shifting, both ef- 
fects occur simultaneously. When the ram pressure ceases, the 
shifted, but bound, ISM reaccretes. 

From Fig.[8]we determine the final post-stripping mass Mgnai 
of the galaxy as the mass of bound ISM at the final simulation 
time of 2 Gyr: Mgnai - 0.71Md4SM- The stripping radius r,,,.,^ 
is trickier to measure, since all the reaccreting material is not 
back in the disk at the final time simulated. We estimate it to be 
r,,„p ~ 6 kpc. Initially in Fig. [8] about 5% of the 12 000 ISM 
particles are, due to relaxation processes and scattering on spiral 
arms, already outside the region r < 16 kpc, |z| < 1 kpc before 
the ram pressure starts to play a role. But all these particles have 
negative total energy. 

The upper row of Fig. |9] shows the ICM density distribution 
in the surroundings of the galaxy passing the densest parts of the 
cluster at times 1.59- 1 .67 Gyr. The particles are displayed as 
filled circles with radii equal to their SPH smoothing lengths /z's, 
where the gray shade corresponds to the density. Due to the su- 
personic motion, we see a bow shock forming in the ICM. It is 
most pronounced when the galaxy passes through the very center 
of the cluster The bottom row of Fig. |9] shows the ISM density. 
Due to large differences in the density between the disk and the 
tail, the grey-scale is logarithmic. Note the growth of ISM par- 
ticles when they are stripped from the disk. In the Appendix, 
Fig.|22]shows the ICM and ISM density distributions for a sim- 
ulation with A^icM - 480 000. 

Figure [10] shows the velocity vectors of the ICM particles 
interacting with the ISM disk. At the surface of the disk, ICM 
decelerates from more than 1000 km s"' to less than 100 km 
s '. We can also see the deflection of the ICM flow to the sides 
along the shock, and a low-velocity ICM in the shadow behind 
the galaxy. In Fig.|23]of the Appendix the ICM velocity field for 
the case of A^icm = 480 000 is shown. 

6.2. LM galaxy in various clusters 

As we have seen in Fig. [8] the stripping of the ISM is a dynami- 
cal process depending strongly on the shape of the ram pressure 
peak. Therefore, we present the results of a set of twenty-five 
simulations of the LM-type galaxy crossing clusters with vary- 
ing distributions of the ICM. 

6.2.1. Results from the simulations 

Multiplying the standard values of the parameters /?c,icm and 
Po.iCM by factors of 0.25, 0.5, 1, 2, and 4, we study the effects 
of a varying width and height of the ICM density peak on the 
stripping results. Fig.[TT]illustrates the ram pressure profiles cor- 
responding to the ICM distribution with varying widths but a 
constant height (left), and varying heights but a constant width 
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Fig. 9. Density of tlie ICM in tlie surroundings of the LM-type galaxy (top) and of the ISM (bottom). A bow shock forms in front of the disk. 
Individual particles are displayed as filled circles with radii equal to their SPH smoothing sizes h's. The scale of the ISM density is logarithmic 
due to large differences in the disk plane and in the tail of stripped particles. The time in Myr is given in the left lower corner of frames. 
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Fig. 10. The flow of ICM particles around the LM-type galaxy. Vectors of the relative velocity of ICM particles interacting with the ISM disk are 
displayed. 



(right). Fig.[T2]shows the stripping results of nine runs from the 
above set of twenty-five simulations. Individual panels are like 
those in Fig. |8]and correspond to po.icivi = 0.25, 1,4x4 10"^' 
cm"'' (from left to right) and .RcICim - 0.25, 1, 4 x 13.4 kpc 
(from top to bottom). 

The main stripping event is rather short in all the cases 
and takes 100 to 200 Myr Going from the upper left corner 
to lower right of Fig. [12] the stripping expands to more and 
more inner parts of the galaxy. It means that both increasing 
width and height of the ICM density peak support the stripping. 
Nevertheless, the process of stripping itself is different in either 
direction: 

a) following a row, i.e. for a fixed value of the ^c.icm, the strip- 
ping of the ISM deepens with growing value of the central den- 
sity. In the upper row, the ram pressure peak is narrow, with 
^c.icivi = 3.4 kpc. If the peak is not high (left), most of the ISM 
remains inside the |z| < 1 kpc layer. About ~ 12% of the ISM has 
been shifted out of the layer, but almost all of it remains gravi- 
tationally bound to the galaxy and later will be reaccreted. The 
wiggles seen on the curves are z-oscillations of the reaccreting 
material around the disk plane. If the ram pressure pulse is high 




R [kpc] R [kpc] 



Fig. 11. Central views of the ram pressure peaks with a fixed po.icM = 
4 10"' cm"' and R,^cM = 0.25, 0.5, 1, 2, 4 x 13.4 kpc (left), and a fixed 
RcicM = 13.4 kpc and po.icM = 0.25, 0.5, 1, 2, 4 x 4 10"' cm"' (right). 



enough (right), it kicks the disk abruptly, which accelerates the 
outer ISM directly to Vg^c and only shifts the inner ISM out of 
the disk plane. Consequently, ~ 52% of the ISM moves out of 
the |z| < 1 kpc layer, but more than a half of it, i.e. ~ 30% of 
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Fig. 12. Stripping of the LM-type galaxy crossing clusters with varying ICM distribution: from left to right: po.icM = 0.25, 1,4x4 10"^ cm"^, 
from top to bottom: Rc.icm = 0.25, 1, 4 x 13.4 kpc. See Fig. [8]for details. Simulations in the same row thus have a constant value of the /?c,icm 
parameter, and simulations in the same column keep the po.icM value. The central panel corresponds to Fig. [8] 



the original amount of the ISM, remains gravitationally bound 
and is later reaccreted. It may be seen as a deep notch with the 
deepest part at T ~ 1.65 Gyr. 

b) Following a column, i.e. for a fixed value of the central den- 
sity, the stripping begins in wider clusters at earlier times, since 
the value of the ram pressure needed to affect the outer ISM 
shifts to greater distances from the cluster center. Thus, the min- 
imum of the ISM mass in the evaluation zone |z| < 1 kpc broad- 
ens. Broad clusters, with Rc,icm - 53.6 kpc, are not only much 
more effective in removing the ISM particles outside the |z| < 1 
kpc layer, but the more extended ram pressure pulse is also able 
to accelerate more ISM and remove it from the parent galaxy. 

Table [3] quantifies the stripping results of the set of twenty- 
five simulations and adds values of parameters characterizing in- 
dividual runs. One sees that in a small (low and narrow) cluster, 
poor in ICM, where po.icM - 1 10"^ cm"^ and /?c,icm = 3.4 kpc, 
the LM-type galaxy is stripped of at most 1 % of its ISM. It hap- 
pens between the orbital time 1.53 and 1.64 Gyr. On the other 
hand, in a large (high and wide) cluster with a lot of ICM, where 
Po.iCM = 1-6 10"^ cm"^ and .Rc.icm = 53.6 kpc, the stripping of 
85% of the ISM is achieved between orbital times 1.37 and 1.56 
Gyr. 

Figure [T3] summarizes the results of the set of simulations. 
The general trends of galaxies crossing wider and higher ICM 
density peaks losing a growing amount of the ISM, as well as 
the reaccretion weakening, are clearly visible. In clusters with 
a constant width, the stripping starts within a short time range, 
while it occurs earlier in wide clusters, where galaxies fall faster 
towards the cluster center and the ICM has a greater extent. 



6.2.2. Comparison with tine criterion of Gunn & Gott 

Focusing in Table |3] on values of the column density of the en- 
countered ICM, we note that the passages through clusters with 
a similar amount of the encountered ICM show a similar amount 
of the stripped ISM. Figure [T4l compares the process of strip- 
ping in five clusters with Eicm ~ (5 - 6) Mo/pc" which are in 
Table [3] marked in boldface. In their case, Ms,rip = 22%, 26%, 
29%, 27%, and 26%. However, the stripping history differs: in 
a narrow cluster with a high ICM density peak, there is a lot of 
ISM shifted out of the \z\ < 1 kpc zone, but this is reaccreted 
later. In a broad cluster with a rather shallow density peak, more 
ISM stays near the galaxy symmetry plane and one sees less of 
the reaccretion. 

The two extreme cases substantially differ in the peak ram 
pressure; in the first case, it is almost 16 x higher compared to 
the second case. Gunn & Gott (1972) criterion (Eq. [TJ would 
predict a completely different amount of stripping in these cases. 
However, as we argue later, the final stripping depends on the 
column density of the encountered ICM, which is similar in the 
two cases. In Fig.[T9]a comparison of the simulation results with 
Gunn & Gott (1972) predictions can be seen. 

In the left panel of Fig. [15] the dependence of the remaining 
gas mass fraction Mpnai = 1 - Ms,np on the extent of the ICM 
distribution in the cluster is plotted. For increasing /?c,icm, the 
degree of stripping becomes less and less dependent on the size 
of the ICM distribution (characterized by the /?c,icm itself) and 
approaches a constant, which depends on the value of pojcM- 
This depends on whether or not the ISM is accelerated during 
the ram pressure pulse to the escape velocity. The critical size 
of the active region AR,. is AR,. = Vesc v SisM/Pra'L' where p"^^^ 
is an effective value of the ram pressure. Here we take p'^!!^^ = 
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Table 3. Summary of stripping results of the set of twenty-five simulations with the LM-type galaxy crossing clusters with varying values of 
the Rcicu and po.icM parameters. Stripping radius r„„p of the ISM disk, mass of the stripped ISM Ms,rip, minimum mass of the ISM within the 
evaluation zone, and mass of the reaccreted ISM Moj.„ are stated. M5,„,, + M„„„ + Maccr = 100%. Values of peak ram pressure p™ , column density 
of the encountered ICM Sicm, and mass of the ICM within central part of clusters Micm are added. 
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Fig. 13. Stripping of the LM-type galaxy in central parts of clusters with: (left pair) fixed Rc.icu = 13.4 kpc and varying Po.icm: 0.25, 0.5, 1, 2, and 
4x4 10"' cm"' (from top to bottom) and (right pair) with fixed po.icM = 4 10"' cm"' and varying Rc,icm'- 0.25, 0.5, 1, 2, and 4 x 13.4 kpc (from 
top to bottom). Evolution of the ISM mass bound to the galaxy or staying within r < 16 kpc and |z| < 1 kpc about the disk plane is displayed. 
Reaccretion of the stripped ISM is observable. 



/7™V2. With SisM = 1021 cm"2, ve^e = 400 km s"i, = 6000 
cm"-^ km^ s"2, and v = 1300 km s"', we get AR^ « 60 kpc. 
This corresponds nicely to Rc.ism - ARc/2 = 30 kpc (see Fig. 
[TSl l. where Monai loses its dependence on R^jsm and converges 
to the Gunn & Gott (1972) prediction given with Eq. [T] Thus 
our results give much less stripping compared to Gunn & Gott 
(1972) in clusters with a narrow ICM peak, while converging to 
its predictions in extended ICM clusters. 

The convergence of the stripping level is apparent from 
Fig. [15] central panel, which depicts the final ISM mass as a 
function of the column density of the encountered ICM. Clearly, 
the higher the Sicm, the more the galaxy is stripped; however, 
the dependence is not linear The decline in the curve corre- 
sponds to the increasing duration of the ICM-ISM interaction 
and to the resulting stripping of more and more strippable ele- 
ments. The final saturation than equals the Gunn & Gott predic- 



tion when all the elements with /,.(,;, < p'"^j^ are stripped. The 
shape of the dependence is unique for the given orbit, galaxy- 
type, and the cluster potential. Then, knowing from observations 
the galaxy's velocity and mass, along with the distribution of the 
cluster DM and Zicm, one can estimate the stripping and find 
limits on the Gunn & Gott criterion. Figure [15] right panel, fur- 
ther shows the reciprocal of the ISM final mass as a function 
of SicM for LM-type galaxy. We try to fit the dependence with 
^ final ~ ^^iCM + 1, where k ~ 0.09. The quantity Mj/^^, in fact 

Con-esponds to M:sMJnil IMisMJinal = MHI,expecled/MHI,ohsen-ed, 

whose logarithm then gives the deficiency DEF introduced by 
Solanes et al. (1996). 
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Fig. 15. Left: Final ISM mass fraction Mflnai as a function of /?c,icm for the LM-type galaxy. The lines correspond to different values of the po.icm: 
1, 2, 4, 8, and 16 x 10"' cm""* (from top down); see legend in the central panel. Center: Mflnai as a function of the column density of the encountered 
ICM, SicMi for the LM-type galaxy. Right: Inverse value of Mfi„ai as a function of Ejcm- 
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Fig. 14. Stripping of the LM-type galaxy in different clusters but with a 
similar amount of the encountered ICM, Ejcm ~ (5 - 6) Mo/pc^ . Curves 
correspond to (Rc,icm, Po.icm) pairs in units (kpc, 10"' cm"'). The ISM 
mass bound to the galaxy (left) or staying within r < 16 kpc and \z\ < 1 
kpc layer about the disk plane (right) are displayed. 

6.3. Various galaxy types 

A comparison of the models Lm, EM, and Em in the standard 
cluster is shown in Fig. [16] The LM-type galaxy has a stripping 
radius of about 6 kpc and loses about 29% of its original ISM. 
For Lm and EM-type galaxies the stripping radii move closer to 
the galaxy center, to about 3.9 kpc and 2.9 kpc, respectively. The 
Em-type galaxy is stripped almost completely. As summarized 
in Table 21 the LM, Lm, EM, and Em galaxy types lose about 
29%, 59%, 69%, and 96% of its original ISM during the face-on 
motion through the cluster with standard ICM distribution. 

We may ask why the amount of the stripping is smaller in 
the case of the Lm-type compared to the EM-type, when the 
maximum restoring force is stronger for the EM galaxy except 
the innermost radii (see Fig.|2l left panel). However, in the range 
(1-4) kpc, where the stripping radii of the two cases reside, 
the maxima of the restoring forces are almost equal. Further, the 
surface density of the ISM is lower, and the escape velocity is 
higher in the EM-type compared to the Lm-type. But, as we see 
later, the combination of a lower escape velocity and a higher 
SisM in the Lm-type does not explain less stripping compared to 
the EM-type. 

As shown in the middle panel of Fig. [21 the maxima of the 
restoring force over z-direction occur substantially closer to the 
disk plane in the case of the Lm galaxy, especially at small galac- 
tic radii. Thus, the rise in the restoring force behind the disk is 
much steeper in the Lm galaxy than in the EM model. In both 
cases, the galaxies move through the narrow standard cluster. 
Consequently, the ISM of the Lm-type galaxy is only shifted 
slightly out of the disk by the short ICM wind pulse. On the con- 
trary, the ISM in the EM-type galaxy is more easily shifted to 
higher z's by the same ram pressure pulse, and is thus acceler- 
ated towards the v^sc more effectively. This difference raises an 
interesting physical possibility: If the peak in the restoring force 



Gal. type 




M,„.,p 




^accr 




(kpc) 


(%) 


(%) 


(%) 


LM 


6.3 


29 


49 


22 


Lm 


3.9 


59 


27 


14 


EM 


2.9 


73 


18 


9 


Em 


0.2 


96 


4 






Table 4. Stripping results for different galaxy types in the standard 
cluster. 



is a long way from the midplane, the ram pressure is able to do 
work on disk gas, building up its kinetic energy, before the gas 
encounters the peak of the restoring force. The stored kinetic en- 
ergy could then be sufficient to get (or helping to do so) the gas 
out of the galaxy, even if the ram pressure alone is insufficient to 
do this. 

The difference in the amount of ICM momentum delivered to 
the ISM also has hydrodynamical causes: the higher density of 
the ISM leads to creation of a larger and stronger bow shock, bet- 
ter protecting the ISM against the incoming ICM. That is more 
effectively deflected to the sides, and a fraction of the ICM pres- 
sure is directed to the sides, thereby decreasing the postshock 
pressure. Then the ISM in Lm-type gets less ICM momenta com- 
pared to EM galaxy. 

Figure [T7] displays five snapshots of the Lm and EM galax- 
ies during their passage through the central parts of the standard 
cluster. The ISM disk of the EM-type galaxy is much thicker 
than that of the Lm-type due to its higher bulge-to-disk mass 
ratio. This leads to another possible explanation of the more ef- 
fective stripping of the EM galaxy, where its thicker disk is more 
vulnerable to the ram pressure. 

7. Simulations with varying SPH sizes of ICIVI 
particles 

To test the influence of the number of ICM particles, we vary 
N\cM by a multiplicative factor a. Then, since the SPH num- 
ber of neighbors is fixed, the smoothing lengths evolve as h' - 
h Q-"'^^. As the number of ICM particles in the simulation in- 
creases, the hydrodynamical interactions are treated in a better 
way, since the mutual sizes of ICM and ISM particles approach. 
We have performed a set of test runs with the LM-type galaxy 
crossing the standard cluster with a number of ICM particles 
varying as 480000, 120000, 24000, and 12 000. Figure [H] dis- 
plays the results - the time evolution of the ISM mass bound to 
the galaxy (top), and the ISM mass located within the r < 16 
kpc, |z| < 1 kpc zone (bottom). Although the number of ICM 
particles differs by a factor 40 in the extreme cases, which cor- 
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Fig. 16. Evolution of the stripping of the Lm, EM, and Em galaxy types (from left to right) during their crossing of the standard cluster. 
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Fig. 17. Edge-on snapshots of the stripping of the Lm (top) and EM (bottom) galaxies. The ISM dislc of the EM galaxy type is clearly thicker in 
comparison with the Lm-type before they enter the ICM distribution. 



responds to about a 3.4-fold difference in the SPH particle sizes, 
only a small difference in the final stripping results is observed 
in Fig. [18] 

Since all particles meeting the condition |r,j| < max(/!,, hj) 
are used to evaluate the pressure gradient at the disk, the ISM 
particles still feel the static pressure of the ICM behind the bow 
shock, even though the SPH sizes of the ICM particles may be 
much greater than those of the ISM. This view is supported 
by Fig. [18] since changing the number of ICM particles does 
not alter the stripping rate. The dominant process governing the 
sweeping in the case of supersonic ISM-ICM interaction is the 
ram pressure. 

Figure[T8]also shows that, due to smoothing of the ICM peak 
by increasing sizes of the particles, shifting of the ISM out of 
the |z| < 1 kpc evaluation zone is steeper. Apart from changing 
sizes of the ICM particles, the growing value of the NycM causes 
different mass ratios of individual ICM vs. ISM particles: 10.3, 
5.2, l,and 0.26. 



8. Analytical model 

After presenting our numerical results, we now consider the in- 
teraction of the galactic ISM with the ICM from an analytic point 
of view. 

8. 7 . The equation of motion of an ISM eiement 

We describe the trajectory of an individual ISM element origi- 
nally sitting in the disk plane, as it moves under the influence of 




,4 1,5 1.6 1.7 1,8 1.9 2 1.4 1.5 1.6 1.7 1.8 1.9 2 
T [Gyr] T [Gyr] 

Fig. 18. Simulation results of the LM-type galaxy crossing the standard 
cluster with a varying number of ICM particles. Total mass of the bound 
ISM (left) and of the ISM located within r < 16 kpc, |z| < 1 kpc zone 
about the disk plane (right) are displayed as functions of the orbital 
time. The vertical line indicates the instant of the passage through the 
cluster center. 



the cluster gravity forces, the restoring force of the galaxy, and 
the ICM pressure force. 

As shown in the simulations, a bow shock forms when the 
incoming ICM arrives at the galaxy's disk. Consequently, the 
ICM decelerates at the shock, transferring part of its momentum 
to the galaxy, and flows along the shock envelope around the 
disk. The disk itself feels the pressure of the bow shock and thus 
the ISM is pushed out of the disk plane, against the galaxy's 
restoring force. 

Disregarding the above hydrodynamical effects and consid- 
ering the face-on case, the equation of motion (EOM) for the 
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Fig. 19. The final mass fraction Mfi„.^\ of the gas that remains in the disk 
after the stripping events as a function of the maximum value of the ram 
pressure p,,,,,,. The symbols (see legend in Fig.|15l> denote the results of 
the SPH calculations for five groups of clusters with increasing central 
density (from left to right), while the gray curves are the corresponding 
predictions obtained by integration of Eq. |6] The thick solid line refers 
to the criterion of Gunn & Gott (1972), Eq.0 



movement of an ISM element in the z-direction, i.e. vertical to 
the galactic plane, is 



; Picm(v - Vour) — T- -^iSM, 

dt oz 



(6) 



where v is the galaxy's orbital velocity. The rightmost term is 
the gravitational restoring force at the current z-position of the 
ISM element, whose current velocity is Vo„t with respect to the 
galaxy's rest frame. The ISM element feels a reduced ram pres- 
sure because of the relative velocity of the galaxy, and the ICM, 
V, is reduced by Vgut- When Vgut exceeds the local value of the 
escape velocity v^scif, z), the element is stripped. 

The solution of this equation, together with dz/dt - Vo„,, 
gives the speed and the trajectory of the ISM element in the 
gravitational field of the galaxy and under the influence of the 
changing ram pressure. From Eq. |6]one notes that, apart from 
the actual value of the ram pressure reduced by the local restor- 
ing force, the history of the encountered ICM also matters. In the 
disk of a galaxy, one applies this EOM to the ISM at each galac- 
tocentric radius; and by numerical integration, one can follow 
the fate of each annulus of the disk until it is stripped. The strip- 
ping radius r,,,.,^ is a radius above which Vo,„(r) exceeds the value 
of the escape velocity in the galactic plane. The final mass is then 
determined as the mass of the ISM enclosed within the stripping 
radius, according to the density profile in the disk (Eq.O. 

Figure [19] compares these predictions of Eq. |6] with the re- 
sults of the SPH simulations, depicting M^nai as a function of the 
peak value pram of the ram pressure. One notes that the curves 
from the simple EOM model match the results obtained by the 
SPH simulations rather well, in particularly for narrow clusters 
with /?c,iCM = 3.4, 6.7, and 13.4 kpc. For large clusters, the 
simple model underestimates the stripping severely. This may 
be because in Eq. |6] we neglect the deformation of the stripped 
ISM disk, which would affect the restoring force, or none of the 
hydrodynamical effects is taken into account. Nonetheless, this 
very simple approach reproduces the dependence on both the ex- 
tend and the density of the cluster quite well. 

Figure [19] also shows that the criterion of Gunn & Gott 
(1972) constitutes the limit for maximum possible stripping. The 
final mass fractions predicted from the simple approach con- 
verge towards the estimates based on Eq. [T] at high values of 
^c.iCM and po.iCM, but they always remain higher This is be- 



cause Gunn & Gott (1972) assume a constant operation of the 
peak value of the ram pressure and compare it with the maxi- 
mum restoring force. Therefore, at all radii of the disk, where 
the (peak) ram pressure exceeds the maximum restoring force, 
all the ISM elements can be accelerated to the escape veloci- 
ties. However, in the integration of the EOM, an ISM element 
may cross the z-distance of the maximum restoring force, but 
the peak of the ram pressure may be over before the element 
reaches escape velocity. 

8.2. The impulse approximation 

Since our principal intention is not to follow the exact trajectory 
of the stripped elements but only to identify them and thus ob- 
tain the stripping radius (or the final mass), let us introduce an 
approximative solution of the EOM. 

If we consider only those regions where the gas is clearly 
removed from the galaxy disk, we can obtain a first-order so- 
lution by considering the balance of the momentum. Neglecting 
the restoring force in these outer regions, the momentum per unit 
surface area attained by the stripped gas of surface density Zism 
and velocity Va/ter is 



after ~ ^after 



2^iSM - I Picmv dt. 



(7) 



This quantity Pafier can be considered as the integrated history 
of the ram pressure: 



af ter 



= 2 



ICM- 



/ Picmv^* 
/ PiCMvdt 



-ICm(v)picm 



(8) 



where Xicm - J PiCM'^dt is the column density of the encoun- 
tered ICM, and {v}p^^^ is the averaged velocity v along the 
galaxy's orbit weighted with the local volume density picM- 
Then Eq. |2]becomes 



Vafter - {v)p 



■'ICM 
-ISM 



(9) 



Of course, the restoring force f^esr can only be neglected in 
the stripped outer parts of the galaxy. A better approximation 
of the EOM is to consider the momentum transferred from the 
inflowing ICM to the ISM element after the element passed the 
position z with maximum value of the restoring force. Thus, one 
regards the acceleration of the ISM only when the ram pressure 
exceeds this maximum. Then, at a given radius. 



fter 



J (PlCMV^ - ^Imax ^ISm) dt 



(10) 



-ISM 



where the integration is taken over a time interval (fi, f2) when 
Pram > fre st,ma.x, i-C- as long as the Gunn & Gott (1972) criterion 
(Eq. [D is fulfilled. Since the actual value of the restoring force 
at the position (r, z) is less compared to its maximum value, this 
obviously overestimates the binding of an element to the galaxy 
plane. 

Figure|20]shows a comparison of the results of the SPH sim- 
ulations and from the approximative solutions by Eqs.l9land[T0l 
which are expected to bracket the actual equation. As before, the 
stripping radius rsirip is the radius above which Vafterif) exceeds 
the value of the escape velocity in the z = plane. Then the final 
mass is determined as the mass of the ISM enclosed within the 
stripping radius, according to the disk's density profile (Eq.[3]l. 
One notes that the curves from Eq. [9] which neglects the restor- 
ing force, indeed overestimate the stripping, while curves based 
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Fig. 20. Numerical results versus analytic predictions of Eqs.[9](left) and [TO] (right) of final remaining gas mass fraction Mflnai after the stripping 
events as functions of the maximum value of the ram pressure /?,„,„. The symbols for individual simulations are the same as in Fig.[T9] The thick 
solid line corresponds to Gunn & Gott (1972), Eq.[T] 



on Eq. \TU\ which consider only the momentum gained as long 
as the ram pressure exceeds the maximum restoring force, un- 
derestimate the stripping. In any case, the departures from SPH 
results are larger than with the predictions of the EOM (Eq. |6]), 
displayed in Fig. [19] 

Of course, comparing the full numerical treatment with the 
predictions of the analytical approach may have several caveats. 
The most important is that all the hydrodynamical processes, 
such as the effect of the bow-shock and size of the galaxy 
shadow, are disregarded in the solution that describes the dynam- 
ics of a single ISM element. Furthermore, in the analytic calcu- 
lations, the density profile of the ISM is assumed to be constant, 
whereas in the full simulations, processes of relaxation and for- 
mation of spiral arms, as well as the progressive stripping, can 
change the surface density profile of the gas disk. Such changes 
consequently affect the restoring force and the escape velocity 
profiles. 

Nonetheless, it is interesting that the simple conservation 
of momentum (Eq. |9]l, which on one hand neglects the actual 
effect of the restoring force and the element's shifting out of 
the disk and on the other overestimates the escape velocity and 
changes of the ISM disk due to stripping, can give a good or- 
der of magnitude estimate of the stripping, when the variation 
in the ICM density along the galaxy's path is to be taken into 
account. This corresponds to the above fit of simulation results 
with M final ~ (k'LfCM + 1)"' : for higher 2/cm, Va/ier grows, turn- 
ing more and more elements to higher than escape velocity v^sc- 
Therefore Mfi„ai decreases with growing E/cm- 

The Gunn & Gott (1972) formula, computed from the max- 
imum ICM density, overestimates the stripping and does not al- 
low for the effect of the finite size of the cluster, i.e. the finite 
time of interaction. However, in cases with high values of /?c,icm 
and po.iCM, the SPH simulations show an even larger stripping 
than the Gunn & Gott (1972) estimate. This may be a con- 
sequence of uneven gas distribution in the galaxy disk where 
spiral-arm inhomogeneities contribute to the ISM removal, as 
also shown by Quilis et al. (2000). 

9. Comparison to other simulations 

9.1. Abadietal. (1999) 

Performing 3D SPH/N-body simulations, Abadi et al. (1999) 
studied the stripping of spiral galaxies by operating a constant 
flow of ICM particles in a simulation box of size 60 kpc x 60 kpc 
X 10 kpc homogeneously filled up with the ICM particles flow- 



ing with periodic boundary conditions. Their model is a spiral 
galaxy exposed to a wind of 1000 km s"', 2000 km s ', or 3000 
km s"' and of density: - 3.37 10~^ cm"^, and p,, = O.lpc- The 
ram pressure then ranges from about 300 to 30 000 cm"^ km^ 
s"^. For the higher value of the ICM density, they estimate the 
stripping radius of a spiral galaxy to 4 kpc, i.e. to about 80 % 
loss of the galaxy's diffuse gas mass. 

9.2. Quilis et al. (2000) 

In their simulations employing a high-resolution 3D Eulerian 
code with a fixed grid based on high-resolution shock-capturing 
method, Quilis et al. (2000) follow the interaction between the 
hot ICM and the cold ISM. The stellar and DM components are 
evolved using a particle-mesh code. They model a luminous spi- 
ral galaxy similar to the Milky Way or Andromeda. The ICM is 
modeled as a uniform medium with temperature Tiqm = 10^ K 
and constant densities - 2.6 lO"-' cm"-' and p,, = O.lpc. The 
galaxy moves through the ICM wind at velocities of 1000 km 
s ' or 2000 km s For the Coma-like ICM density and veloc- 
ity, they estimate the stripping radius as ~ 3 kpc. 

9.3. Vollmeretal. (2001) 

Investigating the role of the ram pressure stripping in the Virgo 
cluster, Vollmer et al. (2001) employ a method of sticky particles 
for modeling the warm neutral clouds of the ISM. The effect of 
the ram pressure is included only analytically as an additional 
acceleration on the clouds located at the windward side of the 
gas distribution. Contrary to other simulations discussed in this 
section, they introduce temporal ram pressure profiles in their 
simulations instead of a constant ICM wind. They use a spi- 
ral galaxy model in a wind creating the ram pressure reaching 
maximum values ranging from 1 000 to 10000 cm"-' km^ s"^. 
Galaxies are considered on slightly elliptical orbits not reaching 
the very center of the cluster Due to a high slope of the yS-profile, 
the ICM density falls steeply with distance R from the cluster 
center. At the smallest pericenter used (R ~ 70 kpc), the ICM 
density exceeds about 3 lO"-' cm"-'. In our present simulations, 
galaxies on completely radial orbits are studied. To achieve sim- 
ilar effects of the ram pressure as Vollmer et al. (2001), we have 
set po.icM = 0.1 X Poic'JIf'^ as our standard value (see Sect. [3). 

Vollmer et al. (2001) show that the ram pressure can lead 
to a temporary increase in the central gas surface density and 
that in some cases a strong reaccretion of the atomic gas occurs 
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Fig. 21. Stripping radius rj,,-,,, versus the peak ram pressure p',™,': a 
comparison with simulations by Abadi et al. (1999) - large X; Quilis et 
al. (2000) - large plus; Vollmer et al. (2001) - large squares; Schulz & 
Struck (2001) - large diamonds; models of Roediger & Hensler (2005) 

- dashed line. Our simulations: -Rc.icm = 3.4 kpc - triangles; 6.7 kpc 

- circles; 13.4 kpc - filled circles; 26.8 kpc - squares; 53.6 kpc - di- 
amonds. Solid line shows the prediction using Gunn & Gott (1972) 
(Eq.[B. 

after the stripping event. They find that the scenario of ram pres- 
sure stripping being responsible for the observed HI deficiency 
is consistent with all HI 21 cm observations in the Virgo cluster. 

The difference between a slightly elliptical orbit with the 
pericenter distance Rp and a completely radial orbit with the 
same maximum value of the ICM density encountered, i.e. 

^o^iCM ~ PiCM^^P^' ^"'^ same /?c,icm, is in the effective size 
of the interacting region. The width of the ram pressure peak 
is larger in the former case than in the latter one. From that it 
follows that in our simulations the ISM-ICM interactions take 
place for a shorter time as compared to Vollmer et al. (2001). 

9.4. Schulz & Struck (2001) 

A model dwarf galaxy in 3D HYDRA (SPH-AP^M) simula- 
tions of Schulz & Struck (2001) is placed in a cubical grid 
with 100 kpc edge. The galaxy is surrounded with 80000 uni- 
formly distributed ICM particles. The temperature of the ICM 
is TicM = 4.6 10^ K, and the density of the wind is ~ 7.3 10"^ 
cm"-*. In their face-on models, the galaxy moves with a velocity 
of 1000 km s ' or 2000 km s Then, the stripping radii mea- 
sured from their Figs. 2 and 14 are ~ 7.2 kpc and ~ 5.2 kpc. 

9.5. Roediger & Hensler (2005) 

With a spiral galaxy models and a wide range of ICM condi- 
tions, Roediger & Hensler (2005) perform a parameter study 
of the ram pressure stripping using 2D Eulerian simulations. 
Except for the gaseous disk, all components are treated 
analytically, providing only the gravitational potential. The 
ICM densities range from 10"^ to 10"^ cm and the result- 
ing ram pressure covers a range from 10 to 10 000 cm"-' km^ s"-. 



9.6. Summary 

The dififerent simulations are listed in Table |5] and a compari- 
son is shown in Fig.|2T] All the models shown, except Schulz & 
Struck (2001), use comparable spiral galaxy models. Abadi et al. 
(1999), Quilis et al. (2000) and Roediger & Hensler (2005) show 



Author 


Method 


ICM flow 


galaxy 


Abadi et al. 


SPH 


constant 


spiral 


Quilis et al. 


3D finite diff. 


constant 


spiral 


Vollmer et al. 


sticky particle 


changing 


spiral 


Schulz et al. 


SPH-AP^M 


constant 


dwarf 


Roediger et al. 


2D finite difi'. 


constant 


spiral 


Jachym et al. 


SPH 


changing 


spiral 



Table 5. The list of simulations. 



high stripping efficiency close to the Gunn & Gott (1972) predic- 
tion, which is due to a constant ICM density adopted. There is 
no difference in the stripping radius and Mfi„ai from 2D and 3D 
hydrodynamical simulations and SPH. It seems that the global 
stripping effect on the galaxy does not depend on the details of 
the hydrodynamical approach. 

The results of Schulz & Struck (2001) show the same strip- 
ping efficiency at lower ram pressures (see Fig. 12111. which is, in 
accordance to Gunn & Gott(1972), due to lower disk restoring 
force in the case of dwarf galaxies. 

Our models with low values of 7?c,icm show less stripping. 
This is mainly due to the effect of changing strength and finite 
duration of the ram pressure peak. For high values of /?c,icm in 
large clusters with long interaction times, our models approach 
the Gunn & Gott (1972) prediction and the results of other mod- 
els. It demonstrates that our SPH approach gives similar results 
to SPH simulations of others and to 2D and 3D Eulerian codes. 

Simulations by Vollmer et al. (2001) use the changing ram 
pressure force acting on sticky particles to show results close 
to Gunn & Gott (1972) prediction. However, with the simi- 
lar dependence of the ram pressure on time, we get less strip- 
ping. Simulations with sticky particles probably overestimate the 
stripping efficiency, since this kind of simulation disregards hy- 
drodynamical effects such as formation of the bow shock in the 
ICM in front of the galaxy. It actually shields the galaxy, thereby 
reducing the efficiency of stripping. 

10. Discussion and conclusions 

The simulations performed in the present work allow us to esti- 
mate the consequences of ram pressure stripping in clusters bet- 
ter If a galaxy like the Milky Way (our type LM) passes in face- 
on orientation through the central parts of a cluster similar to 
the Virgo, it loses about 30% of its interstellar gas. This consti- 
tutes a rather large amount of gas loss. According to Roediger & 
Briiggen (2006), the efficiency of stripping only marginally de- 
pends on the orientation of the galaxy for inclinations lower than 
~ 60°. Therefore, we disregard this aspect in this discussion. 

With this level of stripping, we may estimate the fraction 
of gas supplied to the ICM by stripping. About 2 10^ M© of 
gas is lost by a LM-type galaxy. A much less massive Lm-type 
galaxy is more severely stripped, thus roughly providing the 
same amount of gas. EM and Em-types are subject to an even 
more effective stripping, but their gas disks are less massive. Our 
estimate is that the average contribution per spiral galaxy, late or 
early type, with different orientations in a radial orbit, is 1 - 2 
10^ Mq. This means that the total amount of hot gas in the cen- 
tral region of our standard cluster (6 10'" Mq) can be provided 
by 30 - 60 stripping acts. The amount of stripping from past 
to present is influenced with the decreasing gas content of galax- 
ies, which is an obvious consequence of gas consumption by star 
formation: gas-rich galaxies in the past can provide more gas to 
intergalactic space. In contrast, the growing ICM concentration 
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in the center of the cluster leads to higher stripping efficiency 
now. 

Note that only the contribution from galaxies on radial orbits 
is substantial. Galaxies on circular or elliptical orbits do not pen- 
etrate the dense parts of the ICM and thus are safe from being 
strongly stripped. Depending on the fraction of galaxies in radial 
orbits in the cluster, which unfortunately is not well known, the 
stripping could supply a significant amount of gas to the ICM. 
This may explain the high metallicity of the ICM gas. 

Another consequence of our simulations is that we expect to 
find a significant amount of relatively cold gas in any cluster, 
as debris left over from recent stripping events. This gas is not 
yet mixed with the ICM, but forms large diffuse clouds as tails 
behind the galaxy keeping a fraction of its velocity. The mix- 
ing time should be comparable to the free-fall time, e.g. ~ 1 Gyr. 
The density of gas in tails may be more than 10 times higher than 
the local hot ICM density, hence a tail could cause strong strip- 
ping for any galaxy that happens to cross it at a very different 
velocity. This can then provide a nice explanation of the puz- 
zling strong stripping recently observed in many regions where 
the hot ICM gas is not sufficient to strip. For instance NGC 4522 
is apparently stripped at a large distance from the cluster center 
(Crowl & Kenney 2006), where the ICM density is not sufficient. 
A nice example of a cold HI tail almost 100 kpc long provides 
the galaxy CGCG 97-073 in the Coma cluster This and sim- 
ilar tails may be the cause of ram pressure for other galaxies 
arriving at high speed. The surroundings of Virgo M86 galaxy 
shows a complex structure in X-ray emission (Finoguenov et al. 
2004), like a large plume extending in the northwest direction. 
However, ffie influence of such spatially narrow and small ICM 
peaks acting on a galaxy during a short time is estimated to be 
very slow. Also, ram-pressure stripping appears in small groups, 
where the X-ray gas is not detected (Rasmussen et al. 2006). As 
we have shown, ffie crucial eff'ect of the ram pressure stripping 
is the pressure itself, so ffie hydrodynamical effects play only a 
minor role. Therefore, not only is the hot ICM able to strip the 
galactic ISM, but debris from tidal interactions or previous strip- 
ping do so, as well. 

The main conclusion is that the stripping efficiency depends 
significantly on the duration of the ram pressure pulse, which 
in many cases means that the Gunn & Gott (1972) prediction 
overestimates the amount of stripping. This is nicely shown by 
Boselli & Gavazzi (2006, their Fig. 18), where they compare 
the restoring force in spiral galaxies in Coma, A 1367 and Virgo 
clusters to the maximum ram pressure. It turns out that many 
galaxies with normal HI content should have been stripped if the 
ram pressure had acted long enough. Much more modest strip- 
ping is observed, showing that Gunn & Gott (1972) prediction 
is an overestimate, since it neglects the role of reaccretion of the 
ISM. As encountered in Sect. 16.31 in early-type galaxies with 
high bulge-to-disk mass ratios, the maximum of the gravitational 
restoring force is more displaced from the disk plane than in late 
types. This means that the ram pressure can accelerate the disk 
gas more before it reaches the peak of the restoring force. In prin- 
ciple, this mechanism could allow stripping even when the Gunn 
& Gott condition is not met. The right panel of Fig. [1] shows that, 
especially at larger disk radii, the maxima of the restoring force 
are well off the galaxy midplane. Our simulation results support 
this idea - in Fig. |2T] galaxies crossing clusters with large ICM 
distributions show more stripping than corresponds to the Gunn 
& Gott prediction. 



11. Summary 

Using ffie SPH/tree code GADGET (Springel et al. 2001) 
adapted for interactions of two gaseous phases, the ISM and 
the ICM, we have performed simulations of the ram pressure 
stripping of galaxies on completely radial orbits crossing galaxy 
clusters with realistic profiles of the ICM density. 

We found different amounts of stripping (from complete 
to only marginal) in clusters with various ICM distributions. 
However, when galaxies encounter the same amount of ICM on 
ffieir orbits, the stripping results are found to be the same. Thus, 
for a given orbital type and cluster potential, the column density 
SiCM of the crossed ICM determines the stripping (see Fig.fTSll. 
We have proposed a simple analytic stripping criterion (Eq. |9]) 
taking the Eicm into account. 

Since simulations with different levels of hydrodynamical 
treatment show similar results, we have suggested that the 
ram pressure is the dominant process in sweeping of galaxies. 
Although this also reflects the formation of a bow-shock struc- 
ture in front of the galaxy due to its supersonic motion, we pro- 
posed the important role played by the ram pressure in the inter- 
actions of galaxies with debris structures in clusters. 

We have indicated the limits of the Gunn & Gott (1972) strip- 
ping formula, on one hand overestimating ffie amount of strip- 
ping since completely neglecting the time-dependence of the 
ISM-ICM interaction process, and on the other underestimating 
the stripping by omitting the accumulation of kinetic energy by 
the gas before it encounters the maximum of the restoring force. 
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Fig. 22. Density of the ICM in surroundings of the LM-type galaxy (top) and of the ISM (bottom). N\cm = 480 000. Particles are displayed as 
filled circles with radii equal to their SPH smoothing sizes. The ISM density scale is logarithmic. The time in Myr is given in the left lower corner 
of frames. 
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Fig. 23. The flow of ICM particles around the LM-type galaxy in the case of Mcm = 480 000. 

Appendix 

In Sect. |6?T] Figs. |9] and [TO] show the distribution of the ICM 
density and relative velocity at five locations of the LM galaxy 
when it crosses the standard cluster There, the ICM is repre- 
sented with 120 000 SPH particles. Figuresl22]andl23]depict the 
same situation with Njcm - 480000. Now, ICM particles are 
smaller and the bow shock in front of the galaxy is more pro- 
nounced and ISM in the disk slightly better protected. However, 
as shown in Sect. |7] the stripping results are almost independent 
of the value of Nicm- 



